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ABSTRACT 

We report detection of an extended mid-infrared emission from IRAS 18576+0341 
(AFGL 2298). The object shows a dusty circumstellar shell that has diameter of <; 7" 
at 10.3 and 18.0 (im. The dust nebula shows two emission peaks concentrically elon- 
gated and symmetrically oriented on the opposite sides of the third, central peak, which 
appears to be the central star of the system. The observed mid-infrared morphology 
indicates that the circumstellar dust shell has an equatorially-enhanced material dis- 
tribution, which is a common signature of stellar objects that have experienced mass 
loss. Radiative transfer model calculations suggest that the central star is an extremely 
bright (L* = 10 6 4 -L Q ) star at a distance of about 10 kpc: this object is best described 
as a new luminous blue variable candidate. The circumstellar dust shell seems to have 
been generated by an equatorially-enhanced mass loss process with M > 6.8 x 1O _6 M0 
yr" 1 and M polc /M cq ~ 0.5. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — stars: 
mass loss — stars: individual (IRAS 18576+0341) — stars: individual (AFGL 2298) 
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1. Introduction 

Mid-infrared (mid-IR; 8 — 25/um) imaging provides powerful means to directly probe the distri- 
bution of circumstellar matter, especially circumstellar dust grains. High-resolution (~ 1") mid-IR 
images can reveal the structure at the innermost regions of the circumstellar dust shells (CDSs), be- 
cause the emission is generally optically thin and arises from the warmest (100 — 200 K) dust grains 
in the CDSs. Such mid-IR images of the CDSs have shown that many stars have lost or have been 
losing their mass in an equatorially-enhanced manner, resulting in toroidal or ring shaped CDSs. 
For example, the axisymmetric shaping of planetary nebulae is considered to be initiated by an 
equatorially-enhanced dust-driven wind mass loss at the end of the asymptotic giant branch phase 
(e.g., Ueta, Meixner, Sz Bobrowsky 2000 and references therein). The axisymmetric morphologies 
of the nebulae around luminous blue variables (LBVs) are explained similarly by the presence of 
an equatorially-enhanced CDS (Nota et al. 1995), which may be generated by stellar rotation near 
the Eddington limit (Langer, Garcia-Segura, & Mac Low 1999) or the tidal action of the central 
binary system (Damineli et al. 2000). 

In this paper, we report our discovery of an extended dust shell around an unidentified IR ob- 
ject, IRAS 18576+0341 (IRAS 18576, hereafter), during our observing run to obtain high-resolution 
mid-IR images of CDSs around evolved stars. After reviewing previous observations in §2, we 
present the observed images in §3. The results of the dust radiative transfer model calculations are 
presented in §4, followed by a discussion on the nature of the source in §5. Finally, our conclusion 
is presented in §6. 

2. The Object: IRAS 18576+0341 

IRAS 18576 was first detected as an IR source by the Survey Program of Infrared Celestial 
Experiments sensor and the Far-Infrared Sky Survey Experiment telescope (Price, Murdock, <fe 
Shivanandan 1981). Their successor program, the Air Force Geophysical Laboratory (AFGL) 
Infrared Sky Survey (Price & Murdock 1983), confirmed the detection of this IR source and 
designated the object as AFGL 2298. Infrared Astronomical Satellite (IRAS; IRAS Explanatory 
Supplement 1988) also detected it as a point source and assigned the IRAS ID (IRAS 18576+0341). 
From its IRAS color, IRAS 18576 was classified as a region V object (Van der Veen & Habing 
1988) and was suspected to be a planetary nebula (PN), which is consistent with the source's 
very red nature (class H) revealed by the IRAS Low Resolution Spectra (LRS; Volk & Cohen 
1989; Kwok, Volk, & Bidelman 1997). The LRS also suggested the oxygen-rich (O— rich) nature 
of this object due to its possible silicate absorption feature at 10/xm. Garcfa-Lario et al. (1997) 
obtained near-IR photometric data of the source and classified it as a post-AGB object. However, 
no optical counterpart has been identified (Hrivnak, Volk, & Kwok (2000) report their unpublished 
photometry of V > 23.2). 

Zuckerman & Dyck (1986) were unsuccessful in their search for CO line emission toward 
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IRAS 18576 with the National Radio Astronomical Observatory (NRAO) 12 m telescope. The 
authors, nevertheless, classified the object as an oxygen-rich source based on its IRAS color. No 
detection of 22.2 GHz H2O maser emission was reported in a survey conducted at the Owens Valley 
Radio Observatory (Zuckerman & Lo 1987). Becker et al. (1994) detected 5 GHz radio continuum 
emission in the direction of IRAS 18576 in their Galactic Plane Survey using the Very Large Array 
(VLA) and catalogued the source as GRSR5 37.278—0.226. Based on the positional coincidence, 
they determined the matching probability of the radio source to IRAS 18576 to be 98.4% (White, 
Becker, & Helfand 1991; Becker et al. 1994), suggesting the source to be a PN based on the IRAS 
flux density ratios. In the NRAO VLA Sky Survey (NVSS; Condon, Kaplan, & Terzian 1999), 1.4 
GHz radio emission was measured despite the confused background. 

Hrivnak, Volk, & Kwok (2000) observed IRAS 18576 with the Short Wavelength Spectrometer 
(SWS; de Graauw et al. 1996) on board the Infrared Space Observatory (ISO; Kessler, et al. 1996). 
The spectrum shows an IR excess peaking near 27/xm with unidentified IR (UIR) features at 3.3, 
6.2, 7.7, 8.6, and 11.3 ftm, which are often attributed to polycyclic aromatic hydrocarbons (PAHs; 
Allamandola, Hudgins, & Sandford 1999 and references therein). The ISO detection of carbon- 
rich (C— rich) material, together with the previous O— rich suggestion by the IRAS LRS, makes the 
source both C/O— rich object. The authors provided several possibilities for the nature of IRAS 
18576, which include multiple central sources (H II regions or PNs) and a source with multiple 
shells with differing composition. Table 1 summarizes previous photometric observations of IRAS 
18576. 



3. Observations and Results 

3.1. Mid-IR Imaging 

We obtained images of IRAS 18576 using the University of Arizona/Smithsonian Astrophys- 
ical Observatory Mid-IR Array Camera (MIRAC3, Hoffmann et al. 1998), at the NASA Infrared 
Telescope Facility (IRTF) on 1999 June 7 under photometric conditions. The array is a Boeing 
HF-16 arsenic-doped silicon blocked-impurity-band hybrid array and has a 128 x 128 pixel format. 
The pixel scale was set to 0".33 (42" x 42" field of view), which ensures a Nyquist sampling of the 
diffraction-limited PSF of the telescope. The object was observed with 10% (= AA/A) bandwidth 
at 10.3 and 18.0/um. For flux and point-spread-function (PSF) calibration, we observed a Her (a 
CGS3 standard; Cohen & Davies 1995) before and after the object to check for variations in the 
PSF. Measured PSF size (FWHM) was l".05±0".06 for 10.3/xm and l".53±0".05 for 18.0/xm. With 
an east-west telescope nod throw and a north-south secondary chop throw, we mosaiced a final im- 
age of size about 25" across. Flat-fielding and bad-pixel-masking were applied to sky-subtracted, 
co-added images before mosaicing. Each individual image was subdivided by 4 x 4 during mosaicing 
to make the pixel scale of 0".0825/pix for accurate registration. The total integration times of the 
final images are 720 and 400 sec respectively at 10.3 and 18.0/um, resulting in 1 a rms noise of 14 
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and 94 mJy arcsec -2 . Absolute flux calibration errors are estimated to be approximately 10%. A 
full description of the nod-chop procedure, data reduction, and flux calibration processes can be 
found in Meixner et al. (1999). 

Figure 1 shows the reduced images of IRAS 18576 at 10.3 and 18.0^772. The dust nebula 
appears roughly circular with about 7" diameter and there are local emission peaks in the inner 
region of the nebula. At 10.3^m (Figure la), three emission peaks are clearly seen and are aligned 
to form a straight line whose position angle (east from north) is —3.5°. We shall refer to the peaks 
as northern, central, and southern peaks according to their relative positions. The central peak, 
located near the geometrical center of the nebula, is unresolved and very likely represents the central 
star of the system. The northern and southern peaks are 2" to 2". 5 away from the central peak and 
more extended than the central peak. The arc-like elongation of the northern and southern peaks 
seems to encircle the central peak. At 18.0^m (Figure lb), the morphological characteristics are 
identical to those at 10.3/Lim except that the central peak is less pronounced at 18.0/xm. The total 
specific flux of the nebula is 16.7 Jy and 205 Jy respectively at 10.3 and 18.0^m: the northern peak 
accounts for ~ 70% of the total emission while the rest is shared by the southern peak (~ 25%) 
and the central peak (~ 5%) in both wavebands. Emission from the central peak would impose an 
upper limit for the direct stellar emission, and we estimated the specific flux from the central star 
at each wavelength following the iterative method described in Hawkins et al. (1995). 

Figure 2 shows the temperature and optical depth maps, which are derived from the source 
images by following the method described by Dayal et al. (1998). Temperature (Figure 2a) is almost 
constant over the entire extent of the nebula: 115 — 135 K with an average of 122 K. There is a single 
temperature peak coinciding with the central emission peak and temperature falls off as the distance 
from the central peak increases. This indicates that dust grains in the central peak have higher 
temperature than dust grains in the surrounding nebula (northern and southern peaks). Thus, it 
is reasonable to assume that there is a single energy source located in the central emission peak 
and that the surrounding material is heated by radiation from the central source. This supposition 
agrees with the behavior of the relative strengths of the emission peaks in the nebula. While the 
relative emission strength of the northern and southern peaks remains constant, that of the central 
peak does change at two wavelengths. The observed decrease in the relative emission strength 
of the central peak at the longer wavelength is a natural consequence at the Rayleigh-Jeans tail 
of the blackbody curve with high photospheric/dust temperature, whereas the blackbody curve 
with typical dust grain temperature (<S 200 K) generally peaks at mid-IR regions and would not 
lead to a change in the relative emission strength. The optical depth maps of the dust shell arc 
both similar at 10.3^m (Figure 2b) and 18.0/wn. They indicate that the nebula is optically thin 
(r max = 0.058 and 0.028 respectively at 10.3 and 18.0/um) and the northern and southern emission 
peaks coincide with the regions in which the optical depth is high. This is corroborated by the fact 
that evidence for the central object can be seen both at 10.3 and 18.0[im. Therefore, we interpret 
the mid-IR morphology of IRAS 18576 as limb-brightened edges (northern and southern peaks) 
of an optically-thin edge-on dust torus surrounding the central star, as seen in, for example, HD 
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168625 (IRAS 18184-1623; Robberto & Herbst 1998; Meixner et al. 1999) and ry Car (Polomski et 
al. 1999). We define the inner radius of the dust torus (r^) to be 2". 4, which is the distance from 
the central star to the northern and southern peaks. In Table 2, we summarize quantities of the 
CDS in IRAS 18576 measured and derived from the mid-IR images. 

3.2. Near-IR Photometry 

To check the variability of IRAS 18576, we obtained JHK' photometry of the object with a 
40 inch telescope equipped with the Near-Infrared IMager (NIRIM; Meixner, Young Owl, & Leach 
1999) at Mt. Laguna Observatory 1 on 1999 June 12 and Nov 14 and 2000 August 5 under clear 
skies. For flux calibration, we used ELIAS IR standard (Elias et al. 1992) and UKIRT faint IR 
standard stars (Casali & Hawarden 1992). The standard stars were observed at similar airmasses 
as the source or at low (~ 1) and high (J> 2) airmasses to perform an airmass correction. Data 
were taken by shifting the telescope with a 9-point dithering pattern. Each exposure was flat- 
fielded to eliminate large pixel-to-pixel sensitivity variations in the detector array and was then 
sky-subtracted before being co-added into a single frame. Flats were created for each waveband by 
exposing the twilight sky. Sky emission maps were constructed by taking a median of all dithered 
frames after unusually low and high pixels had been masked out. Our reduction method generally 
follows the method described in McLean & Teplitz (1996). Table 3 shows JHK' photometric data 
for IRAS 18576, in which all magnitudes have been converted to the ELIAS (CIT) magnitudes 
(Elias et al. 1992) for comparison. 

It appears that near-IR brightness of IRAS 18576 is in decline in the past decade with the most 
recent observation (2000 August 5) marking the dimmest ever. Because we did not detect a similar 
brightness variation in other objects observed on the same nights during the NIRIM observations, 
we eliminate the possibility of any systematic and/or instrumental effects. Therefore, it is likely 
that IRAS 18576 varies its near-IR brightnesses by a few tenth of a magnitude over a couple of 
months and about 0.5 magnitudes over a decade. 

4. Radiative Transfer Modeling 

4.1. 2-D Dust Radiative Transfer 

To model IR continuum emission and mid-IR morphologies, we have performed model calcu- 
lations for an axisymmetric CDS by using a code which solves the radiative transfer equation in a 
fully two-dimensional grid. This code treats dust absorption, reemission, and (isotropic) scattering 



1 Mt. Laguna Observatory is jointly operated by San Diego State University and University of Illinois at Urbana- 
Champaign. 
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in an equatorially-enhanced dust distribution (i.e., a presumed radial and latitudinal density vari- 
ation) with an arbitrary optical depth for a given set of dust species of a single grain size (Meixner 
et al. 1997; Skinner et al. 1997 for more description of the code). The equatorially-enhanced dust 
distribution in the model CDS assumes a two-phased mass loss process, in which the central star 
changes its mode of mass loss from spherically symmetric to axially symmetric with a preferential 
mass loss towards equatorial directions near the end of the entire mass loss phase. As a result, the 
pole-to-equator density ratio at the inner shell boundary can be less than unity. The code requires 
input parameters that can be categorized into three types: stellar, shell, and dust parameters (Ta- 
ble 3). Stellar parameters (the effective temperature, T e g, radius, R*, and distance) determine 
the available energy flux for dust heating. Shell parameters (inner and outer shell radii, ra, min 
and 7"d, max) optical depth along the polar direction at some reference wavelength, t\, and other 
geometrical parameters) specify the dust shell (both physically and morphologically) for radiative 
transfer calculations. Dust parameters (grain radius and absorption and scattering coefficients, 
Qabs and Q S ca) determine opacity for the shell material. 

The use of this dust radiation transfer code is appropriate because direct stellar radiation 
is estimated to be the dominant heating source for the circumstellar dust grains. The total far- 
infrared flux (observed specific flux integrated beyond 20fim) of IRAS 18576 is ~ 3.9 x 10~ n W 2 
m~ 2 . However, we estimate the amount of far- infrared flux expected from dust grains heated by 
Ly a flux (i.e., line emission) to be ~ 8.7 x 10~ 13 W 2 m -2 , following the the method described by 
Zijlstra et al. (1989). Thus, we can safely neglect dust heating due to line emission. 

When available, stellar and dust parameters are taken from previous observations in the lit- 
erature. Otherwise, we determine these parameters through iterative model calculations seeking 
the best-fit to the spectral energy distribution (SED). Once we achieve the best-fit to the SED, 
we further iterate calculations with different sets of shell parameters and inclination angle 2 (#i nc i) 
seeking the best-fit to the SED as well as the mid-IR morphology, which tends to be strongly in- 
fluenced by the geometric parameters. For example, limb-brightened peaks of a dust torus would 
not be resolved with a too high inner shell optical depth and a single shell structure can yield 
different morphologies depending on the inclination angle. The overall shape of the nebula would 
be ring-like if 6* mc i ~ and would be flattened in the direction of the toroidal axis if 6>i nc i ~ 90. 
The final best-fit model was achieved, following the guidelines explained above, after ~ 200 model 
calculations shifting carefully through the parameter space. Because our code does not handle dust 
stratifications at the moment, the final best-fit model was constructed by combining separate best- 
fit models for each of the O— rich and C— rich components of the dust shell, assuming the cavity 
inside the axisymmetric O— rich shell is filled with spherically symmetric C— rich material. 



2 The inclination angle is the sharp angle between the line of sight and the polar axis. 
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4.2. Modeling Results 



Figure 3 shows the best-fit SED of our model calculations. The best-fit SED (thin solid line) 
consists of possibly stratified shells of differing composition: the outer O— rich shell (thin dashed 
line) and the inner C— rich shell (thin dash-dotted line). For comparison, the ISO spectrum (thick 
solid line), the IRAS LRS 3 (thick dashed line), and other photometric data (symbols) are shown. 
Table 3 lists input and derived quantities for the best-fit model. 

We also show the best-fit synthesized images at 10.3 and 18.0/zm (Figure 4). The final syn- 
thesized images were constructed by adding the observed standard star images to the O— rich shell 
model images because the code does not produce images with a resolved central star. The model 
images were convolved with a Gaussian profile whose F WHM is equal to the observed FWHM at 
each waveband before the central star was added. The two emission peaks arise from equatorially- 
enhanced (the pole-to-equator density ratio = 0.5) O— rich dust distribution viewed with 45° incli- 
nation angle. The resulting morphologies are very sensitive to the inclination angle when the shell 
is optically thin and would completely be altered if the inclination angle is changed by ~ 10%. The 
axis of the dust torus is rotated in the plane of the sky at a position angle of 86.5° East of North 
and the near side of the dust torus is to the east from the center. The slight asymmetry seen in 
the shape of each emission peak is due to self-extinction arising from the polar density variation. 



Because we have no prior knowledge about the central star, we must iteratively derive T e g 
and i?* from radiative transfer calculations. Consider an optically thin, isothermal (at Tj) dust 
shell. When dust grains of a uniform grain size are in thermal equilibrium with the ambient stellar 
radiation field, we can relate the luminosity of the star (L*), shell radius (r^), and dust temperature 
by equating the rate of energy absorption and reemission as 



Here, we have assumed a power-law dust grain emissivity (Qabs ^ n ) an d the constant of propor- 
tionality is solely determined by the dust grain properties. Because both L* and r<i are observable 
quantities modulo distance, we can eliminate the distance dependence in eq.(l) and rewrite it as 



5. 



Discussions 



5.1. 



The Central Star and its Distance 




(1) 




(2) 



3 IRAS LRS spectrum does not have an absolute flux calibration and is shown primarily to show the shape. The 
best-fit was determined by comparing model fluxes with the ISO spectrum and other photometric data. 
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in which F* is the flux of the source and r' d is the angular radius of the dust shell in arcsec. The 
dust temperature effectively defines the location of the IR excess peak in the SED. Therefore, eq.(2) 
indicates that we can fix F* by fitting the IR peak location (i.e., Td) by iterating on and dust 
species in radiative transfer calculations because r' d is observationally determined. It is also possible 
to constrain F* by using the observed IR flux of the source (F° bs ). Photometric data obtained by 
ISO (Hrivnak, Volk, & Kwok 2000) and other near-IR observations (Garcfa-Lario et al. 1997) 
yielded F° bs ~ 2.3 x 10~ 7 erg s~ x cm~ 2 . The near-IR part of the ISO spectrum (Hrivnak, Volk, 
& Kwok 2000) suggests the presence of PAHs. However, solely C— rich models yielded a warmer 
dust temperature than can be produced by the detached dust shell of 2". 4 angular radius due to 
generally higher opacity of carbonaceous material. Models with a mixture of amorphous carbon 
and astronomical silicates produced too high dust temperature as well. Our iterative calculations 
concluded that F, ~4x F° bs with astronomical silicates would produce the best-fit mid-IR SED 
to observations. This means that the O— rich dust shell processes about 25% of the stellar flux, 
and it is consistent with our finding that the dust shell of IRAS 18576 is optically thin. 

To recover T e g and R* from the derived F*, we need to determine the distance to IRAS 18576. 
Previous CO observations do not provide any kinematic information of the source (Zuckerman & 
Dyck 1986) probably because of the confused background reported by Condon, Kaplan, & Terzian 
(1999). Recent CO observations failed to obtain any reliable kinematic data due to the lack of clear 
background in the crowded radio field around IRAS 18576 (D. Fong 1999, priv. comm.). Therefore, 
we estimated the distance to IRAS 18576 by determining the total visual extinction (Ay) that is 
required to make our model near-IR SED fit the observed one. The near-IR colors of IRAS 18576 
are very red. Because the dust shell is assumed to process only about 25% of the stellar radiation, 
it is reasonable to explain the red near-IR colors by a severe interstellar extinction. Our best fit 
model indicated that Ay would have to be about 28 to account for the extinction assuming a 
typical total-to-selective extinction ratio (Ry) of 3.1. Such a high Ay value may not be surprising 
because of the Galactic coordinates of IRAS 18576, (I, b) = (37.3°, —0.3°): the line of sight lies in 
the Galactic plane grazing the outskirts of the bulge. According to the full-sky dust maps compiled 
from COBE / DIRBE and IRAS maps (Schlegel, Finkbeiner, & Davis 1998), the Galactic selective 
extinction, E(B — V), along the line of sight to IRAS 18576 is estimated to be 18.25, which equals to 
the Galactic total visual extinction of 56.58. If we assume the Parenago type exponential distance 
scaling of Ay (Parenago 1940), Ay ~ 28 can be achieved if IRAS 18576 is located around 10 kpc 
away. This estimate is uncertain because the assumed distance scaling law of Ay does not account 
for local clumping of reddening material and may locally have a drastically different form. 

Alternatively, the distance to IRAS 18576 can be estimated by iteratively searching for a 
high enough to heat the optically thin dust shell of IRAS 18576 (Figure 2b) to the observed 
temperature (Td = 115 — 135 K) because the distance directly scales the luminosity of the object. 
The optical depth maps have yielded V6ry low T V£llu.6S (T ma x — 0.058 and 0.028 at 10.3 and 
18.0/mt, respectively). With these optical depth, IRAS 18576 would have to be located about 25 
kpc, which yields an unrealistically high luminosity for a star and brings the object close to the 
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edge of the Galaxy in the direction of the object. Thus, we kept increasing the optical depth of the 
dust shell until the toroidal structure would be unresolved while iteratively adjusting the distance. 
After further iteration, we concluded that 10 ± 3 kpc would scale L* high enough so that model 
calculations would yield the dust peak in the observed SED, using the highest optical depth possible 
to reproduce the resolved dust toroidal structure in our models (Table 3 and Figure 4) . Considering 
the above two independent means to estimate the distance, we adopted 10 kpc to be the distance 
to IRAS 18576 because there is no alternate source of data. With the adopted distance of 10 kpc, 
L* is 1.1 x 10 40 ergs s^ 1 (= 1O 6 - 4 L ). 

We can then constrain T e g from the position of the redward slope of the stellar emission peak 
with respect to the observed SED. A lower T c g would shift the stellar peak to the red and would 
effectively increase the stellar emission in the wavelengths of our interest while leaving the dust 
peak virtually unchanged. Thus, the lower limit for T e g- can be obtained because any lower T c g 
than the limit would yield too much stellar flux than observed. We used the 10.3/xm emission in 
constraining T e g because this wavelength range seemed to retain the most pristine stellar emission 
characteristic free from extinction and/or dust feature emission. Further iteration yielded the lower 
T e Q limit of ~ 9000 K, which corresponds to the upper i?* limit of ~ 7OOi?0. However, it is more 
difficult to constrain the upper T e g limit defining the possible range for T e g because the emission 
properties in the optical to near-IR depend more heavily on the distance to the object (i.e., Ay) 
and the location of the C— rich shell that is considered to be situated inside the O— rich dust shell 
(see below). Nevertheless, iterative calculations suggested that T c g ~ 15000 ± 6000 K would be 
a likely range of the effective temperature. The model calculations with this T e g range requires 
Ay = 24 ~ 32 to be consistent with the observed data. 



5.2. Dust Shell Composition and Structure 

The presence of silicates in IRAS 18576 was originally suspected from the apparent absorption 
feature at lOjum (Zuckerman & Dyck 1986; Volk & Cohen 1989) seen in the LRS. This 10/xm 
dip, however, is a false depression caused by the strong neighboring PAH features at 8.6 and 11.3 
fim, because the CDS is observed to be optically thin enough that the central star is visible even 
at the mid-IR (Figure 1). This makes it difficult to identify the silicate species from the shape 
of the 9 — llfim silicate features (Speck 1998 and references therein). The featureless far-IR 
portion of the spectrum also makes it difficult to determine if there exists a particular crystalline 
silicate species. Therefore, we used optical constants for astronomical silicates (Laor &; Draine 
1993) which empirically simulate a mixture of different types of silicates. With iteratively obtained 
~ 4 x F° hs , the model SED fits the overall shape of the ISO SED very well at the long wavelength 
regions (;> 18 fim) and at the 10/Ltm trough (Figure 3, thin dashed lines). 

The model SED arising only from the O— rich shell shows emission "deficits" in the near-IR 
(<; 10/xm) and at the blueward shoulder of the IR peak (10 — 18/xm). To account for the presence of 
the observed C— rich material, we considered an inner C— rich (PAHs) region, which is surrounded 
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by the O— rich shell of angular radius 2". 4. C— rich material is expected to be closer to the central 
star because of the presumed ionization state of PAHs: the comparable strengths of 7.7 and 11.2 
fim emission peaks suggest that PAHs are marginally ionized (Molster et al. 1996). The inner 
C— rich shell adequately compensates the emission deficits (Figure 3, thin dot-dashed lines). The 
C— rich shell does not reproduce each emission line of the UIR features because we have used 
optical constants of amorphous carbon grains of 0.002/im to substitute for the yet unavailable 
optical constants of PAHs. This dust grain size is comparable to the major PAH cluster size of 
300 — 400 C— atoms (Omont 1986) whose presence is typically observed through 5 — 10/im and 
11 — 15/xm emission plateaus (Beintema et al. 1996). Due to lack of spatial information in the 
near-IR, we determined the inner C— rich shell radius through model iterations. 

The deficit at the longer wavelengths (15 — 20[im), on the other hand, seems to be excessive 
because very little PAH emission is expected in this wavelength region. To boost IR emission from 
15 to 20/im, we included amorphous aluminum oxide (AI2O3), whose Q a bs has a local maximum 
around 13/zm, according to the mass abundance ratio consistent with the cosmic abundance (~ 96% 
amorphous silicates and ~ 4% aluminum oxides). The inclusion of aluminum oxide grains is not 
entirely ad hoc; they are considered to be the nucleation cores in the heterogeneous dust nucleation 
theory (Speck 1998 and references therein) and the ISO spectrum shows a number of possible 
narrow peaks in the 11 — 13fim plateau (especially 11.9 and 12.8 /ira), which may be due to 
aluminum oxides (Begemann et al. 1997). Although we did not attempt model calculations with 
any specific silicate species, the presence of an emission plateau near 34^m, which our model SED 
fails to reproduce (due to the use of amorphous silicates), suggests the presence of crystalline 
olivines (forsterites), which have an emission peak in that region. 

Our best-fit model agrees with the general behavior of the observed SED with the 10 /im trough 
adequately reproduced. This method obviously does not conserve flux at the interface between the 
two shells in the sense that the O— rich shell receives stellar flux which should have been reduced 
by the presence of the C— rich shell. However, the very optically thin C— rich shell is estimated 
to require only about 0.5% of the total flux of the central star. Therefore, considering that the 
O— rich shell itself processes only ~ 25% of the total stellar flux, we conclude that our results would 
still hold even if calculations were done in a strictly flux conserving manner. 

5.3. Nature of IRAS 18576+0341 

The derived stellar parameters put the star above the empirical luminosity limit (Humphreys- 
Davidson limit; Humphreys & Davidson 1979) in the Hertzsprung-Russell diagram. IRAS 18576 
would then likely be a luminous blue variable (LBV) or a red super giant (RSG). The derived 
luminosity of IRAS 18576 is unusually high (L* = 10 6 4 L Q ) for a RSG, however, there exist LBVs 
whose L* is comparable to that of IRAS 18576 (e.g., r\ Car and AG Car, Humphreys & Davidson 
1994 and references therein). Although the expected temperature range for the object is rather low 
(T cff = 15000 ± 6000 K) for a typical LBV, it is still higher than that of the coolest LBVs. The 
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fact that T c fj would not be lower than 9000 K seems to favor the LBV interpretation. One of the 
coolest LBV candidates, HD 168625, was found to have T e g in the range of 12000 — 15000 K (Nota 
et al. 1996). Interestingly, this cool LBV is also found to be associated with an extended dust shell 
with the equatorial enhancement (Robberto & Herbst 1998; Meixner et al. 1999). 

There are a number of observational similarities between IRAS 18576 and other LBVs. The 
mid-IR morphology of LBVs is characterized especially by their signature limb-brightened lobes 
with a separate central source (e.g., r] Car, Polomski et al. 1999; AG Car, Trams, Waters, Sz Voors 
1996; HD 168625, Robberto & Herbst 1998; Meixner et al. 1999), which are very similar to what we 
have seen in IRAS 18576 (Figure 1). Morphologically, it is very likely that the central emission peak 
represents the one and only energy source in the system and illuminates the surrounding lobes. The 
limb-brightened lobes strongly suggest the presence of an equatorially enhanced dust distribution 
in the CDS, and it is consistent with what our model calculations assume. This axisymmetric 
geometry has been thought to be a common signature of the LBV nebulae caused by a sharp pole- 
to-equator density contrast (Nota et al. 1995). Although the origin of the equatorial enhancement 
is not yet fully understood, a number of possible mechanisms have been suggested. For example, 
Langer, Garcia-Segura, & Mac Low (1999) considered rotating stars near their Eddington limit 
and showed hydro dynamically that highly non-spherical mass loss would occur. On the other hand, 
Damineli et al. (2000) have spectroscopically confirmed the binarity of r/ Car, and it has also been 
hydrodynamically shown that such binary systems can generate an equatorially-enhanced material 
distribution around the central system in the context of the PN formation (Mastrodemos & Morris 
1999). Whatever the true scenario may be, the observed limb-brightened mid-IR lobes clearly 
showed the non-sphericity of the system, and any dynamical models should be able to reproduce 
such geometry. 

Spectral characteristics of IRAS 18576 are also similar to those of other LBVs. Large IR 
excesses due to dust grains have been observed in LBVs, and the observed mid-IR flux levels of 
IRAS 18576 are comparable to other LBVs: it is approximately a few times brighter than typical 
LBVs (Lamers et al. 1996) but a couple of orders of magnitudes dimmer than t] Car (Morris et 
al. 1999). Both amorphous and crystalline silicates have been identified in the LBV dust shells 
(Waters et al. 1997). The presence of C— rich material in a dominantly O— rich shell has already 
been seen in other LBVs (e.g., Voors et al. 1999). Moreover, IRAS 18576 seems to be varying its 
brightness with the timescale (months to years) and the widths of magnitude (about 0.1 — 0.5) 
consistent with what have been observed in LBVs. 

The observed and derived physical quantities of the dust shell are also consistent with LBV 
characteristics. The maximum extent of the observed mid-IR dust nebula (<; 5a detection) is about 
10", which equals to ~ 0.5 pc at 10 kpc. Considering the fact that the mid-IR observations probe 
only the innermost regions of the dust shell, this value is reasonable for the inner radius of a typical 
LBV nebula (e.g., Nota et al. 1995; Chu, Weis, & Garnett 1999). The total dust mass in the IRAS 
18576 dust torus is estimated to be ~ 0.1M Q . This estimate depends on the outer radius of the 
dust torus, which is assumed to be 6 times larger than the inner radius. The estimated dust mass 
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is nearly a factor of 10 higher than typical LBV dust mass quantities (Nota et al. 1995) but is 
comparable to the estimated dust mass of a similar temperature (100 ± 10 K) component in rj Car 
(Morris et al. 1999). The total mass lost by the star and the rate of mass loss are estimated to be 
about 75M Q and 6.8 x 1O" 6 M yr 1 using the canonical gas-to-dust ratio of 100 and an assumed 
expansion velocity of 50 km s^ 1 (e.g., Nota et al. 1995). Both of these values depend on the outer 
radius of the entire mass loss shell, which is not constrained by observations and is assumed to be 
50 times as the inner shell radius. 

Based on the mid-IR morphology and IR photometric and spectroscopic characteristics, and 
results from radiative transfer calculations, IRAS 18576 is likely a very luminous star surrounded by 
an optically thin, possibly stratified C/O— rich dust shell located at about 10 kpc in the direction 
in which very severe interstellar extinction can occur. Physical parameters of the IRAS 18576 
system derived from radiative transfer calculations suggest that IRAS 18576 is probably a LBV 
candidate. This conclusion derived from the fact that the derived luminosity for IRAS 18576 was 
so high and that other evidence seemed to be consistent with our conclusion. However, given the 
uncertainties involved in the distance determination, it is also possible that IRAS 18576 is a RSG. 
Therefore, it is necessary to observationally uncover the physical characteristics of the central star 
to determine the exact evolutionary status of IRAS 18576. Because IRAS 18576 is not visible in 
the optical due to heavy extinction, we need near-IR spectroscopic information to constrain the 
physical parameters of the central star. We also need to continue obtaining near-IR photometry to 
better characterize the variability of the source. We are planning such follow-up observations. 

Regardless of the exact evolutionary status, IRAS 18576 seems highly likely to be an evolved, 
massive post-main-sequence star. The axisymmetry of the circumstellar shells of evolved massive 
stars has long been indirectly suggested from the strongly non-spherical shapes of optical nebulosi- 
ties associated with LBVs (e.g., Nota et al. 1995 and references therein). Our discovery of the dust 
shell around IRAS 18576 is one of the few direct evidence for the equatorial density enhancement 
in the circumstellar environment of massive post-main-sequence stars (e.g., r] Car, Polomski et al. 
1999; AG Car, Trams, Waters, & Voors 1996; HD 168625, Robberto & Herbst 1998; Meixner et al. 
1999). For low to intermediate mass stars, it is well established that the mode of mass loss is highly 
axisymmetric (e.g., Ueta, Meixner, & Bobrowsky 2000). It now appears that an axisymmetric 
mass loss also operates for massive stars, and therefore, stellar wind scenarios for massive stars 
would be able to produce the observed pole-to-equator density contrast and generate axisymmetric 
CDSs. 

6. Conclusions 

We have discovered an extended mid-IR emission around IRAS 18576 and presented the mid- 
IR images at 10.3 and 18.0 /j,m. The two emission peaks, which are concentrically elongated and 
symmetrically oriented on the opposite sides of the central peak, are typically the results from an 
optically thin, limb-brightened edges of the dust torus viewed rather edge-on. We also present 
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new near-IR photometry of the object which demonstrates its variability. Although the source 
lies in the region of the sky in which a severe Galactic extinction has been observed, we have 
constructed a model using all of the available information. The derived physical quantities for 
the central star and dust shell and variability of the object suggest that IRAS 18576 is a possible 
LBV candidate with L* = 1O 6 - 4 L and T cff = 15000 ± 9000 K. The uncertainties involved in the 
distance determination leave possibilities that the object could be a RSG, and further observational 
work is needed to determine the exact evolutionary state of the star. While the main dust shell is 
composed of silicate dust grains, an additional C— rich (composed of PAHs) shell likely exists inside 
of the O— rich shell. The axisymmetric nature of the dust toroid is characterized by the calculated 
pole-to-equator density ratio of 0.5. While the origin of the equatorial enhancement is still open 
for further discussion, the axisymmetry seems to be a common feature in mass loss from massive 
stars as in the case for low to intermediate stars. 
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Fig. 1. — Observed grayscale images of IRAS 18576+0341 at 10.3/xm (a; top) and 18.0/xm (b; 
bottom) with north being up and east being to the left. The tick marks show relative offsets in 
arcseconds. Contours are spaced by 10% of the peak intensity and their colors are inverted for the 
ease of viewing. The lowest contour is equivalent of 4.3a and 6.7a level of emission respectively at 
10.3 and 18.0 \im. Local peak intensities are, from northern to southern peak, 0.57, 0.49, and 0.38 
Jy arcsec" 2 at 10.3 \im and 6.3, 3.9, and 4.4 Jy arcsec -2 at 18.0 fim. The PSF size (FWHM) is 
indicated by a black circle in each frame at the lower left: l".05 at 10.3^m and l".53 at 18.0/xm. 

Fig. 2. — Temperature map (a; top) and optical depth map at 10.3/^m (b; bottom) of IRAS 18576 
following the display convention of Figure 1. (a) Temperature contours go from 135 to 115 K 
with an interval of 2.5 K. The spurious peaks around the nebula are artifacts occurred during the 
derivation, (b) Optical depth contour interval is 10% of the peak value, 0.058. 

Fig. 3. — The spectral energy distribution of the best-fit model (9000 K, 700-R*, 10 kpc): thin solid, 
dashed, and dot-dashed lines respectively show the total SED, contribution from the O— rich shell, 
and contribution from the C— rich shell with ISM extinction (bottom lines) and without extinction 
(top lines). Observational data are indicated by thick lines (solid — ISO ; Hrivnak et al. 2000, 
dashed — IRAS LRS; Volk and Cohen 1989) and symbols (circles — IRAS photometry; IRAS 
Explanatory Suppliment 1988, asterisks — AFGL photometry; Price & Murdock 1983, crosses — 
near-IR photometry; Garcfa-Lario et al.1997, stars — MSX photometry; Egan et al.1999, squares; 
this observation). 

Fig. A.— Model Images of IRAS 18576+0341 at 10.3/xm (a; top) and 18.0/xm (b; bottom) following 
the display convention of Figure 1. Contours are spaced by 10% of the peak intensity. Peak 
intensities are 0.40 Jy arcsec -2 (at 10.3 fim) and 8.3 Jy arcsec~ 2 (at 18.0 fim). The central 
star (the observed PSF) was scaled to reflect the observed northern-to-central peak ratio at each 
wavelength. 
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Table 1. Previously Obtained Photometry on IRAS 18576+0341 



RA (2000) 


DEC (2000) 




12.0fmi h 


19.8^m a 


25.0fmi h 


27.4^tm a 


60.0/im b 


100.0Atm b 5 GHz c 


1.4 GHz d 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) (mJy) 


(mJy) 


19:00:11.2 


+03:45:46 


81.4 


58.48 


356.9 


424.2 


649.4 


274.5 


< 1661 78.1 


8.1 ±0.7 



References. — a: AFGL survey (Price & Murdock 1983), b: IRAS (IRAS Explanatory Supplement 1988), c: GPS survey 
(Becker ct al. 1994), d: NVSS survey (Condon, Kaplan, & Terzian 1999) 

Note. — See Table 4 for the near-infrared photometric data and Figure 3 for IRAS LRS and ISO spectra. IRAS lOO.Ofim 
data is an upper limit. 
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Table 2. Measured and Derived Quantities of the Dust Shell around IRAS 18576+0341 



Quantity 


10.3 fim 


18.0 fim 


Flux Density (Jy) 


16.7 


205 


Diameter (FWHM) 


6".2 


6".9 


Toroidal Shell Radius 


2" 


.4 


Dust Temperature 


115 ~ 


135 K 


Optical Depth 


& 0.058 


£ 0.028 



-20- 



Table 3. Input and Derived Model Quantities 



Stellar Parameters 




10 6 - 4 L Q 




^cff 


15000 ± 6000 K 




Ft* 


120 - 7OOR 




d 


10 ± 3 kpc 




ISM A v 


24 - 32 




Dust Shell Parameters 




O-rich Shell 


C-rich Shell 


I'd, min (PC) 


0.12 


0.0038 


»*d, max (PC) 


5.8 


0.12 


^incl 


45 ±5° 




T d at r dj min (K) 


112 


512 


^lOSfim, eq 


0.54 


0.036 


^~10.3fj,m, pole 


0.15 


0.036 


Ppole/Peq 3-t 7^ m j n 


0.5 


1 


Mdust (M@/yi) 


6.8 xl0~ 6 


5.7 xl0~ 9 


^exp a 


50 km s" 1 


50 km s- 1 


Tdyn (yr) b 


1.1 x 10 5 


2 x 10 3 


a 


0.01 fmi 


0.002 urn 


Composition 


96% Amorphous Silicates 100% Amorphous Carbon 




4% Alminum Oxides 





a Typical u e xp for LBVs (e.g., Nota et al. 1995). 

k^dyn — r d/ v cxp 

c Abundance % by mass; References: amorphous silicates; Laor & Drainc 

(1993) , aluminum oxides; Begemann et al. (1997); Harman, Ninomiya, & Adachi 

(1994) , amorphous carbon; Zubko ct al. (1996). 
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Table 4. JHK' Photometry of IRAS 18576+0341 



Date 


J 


H 


K' 


Ref. 


1989 May-Jun 


12.0 ±0.3 


8.99±0.03 


6.91±0.02 


1 


1999 Jun 12 


12.21±0.12 


8.94±0.06 


7.07±0.05 


2 


1999 Nov 14 


12.44±0.07 


9.32±0.05 


7.60±0.08 


2 


2000 Aug 5 


12.61±0.51 


10.44±0.16 


8.58±0.11 


2 



References. — 1. Garcfa-Lario et al. (1997), 2. NIRIM observa- 
tions. 
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